% Crisis - Taylor Rule response:
[ rf_taylor, rh_taylor, n_taylor, aggc_taylor, m_taylor, pie_taylor, y_taylor  ] = key_taylor(  );

% Crisis - LSAP response (without adjusting Rn):
[ Rf_lsap, rh_lsap, n_lsap, aggc_lsap, m_lsap, pie_lsap, y_lsap ] = key_lsap(  );

% Figure
%----------------------------------------------------------------------------
%List of Figure 
%----------------------------------------------------------------------------
set(groot, 'DefaultLegendInterpreter', 'latex')

figure('Name', 'Federal Funds Rate and Real Borrowing Rate');
plot(truc, rf_taylor, truc, rh_taylor, truc, Rf_lsap, truc, rh_lsap);
leg = legend('$R^f$ - Taylor', '$r^h$ - Taylor', '$R^f$ - LSAP', '$r^h$ - LSAP', 'location', 'best'); 

figure('Name', 'Real balance of Reserves');
plot(spline_truc, n_lsap, spline_truc, n_taylor);
leg = legend('$n$ - LSAP', '$n$ - Taylor', 'location', 'best');

figure('Name', 'Aggregate consumption');
plot(spline_truc, aggc_lsap, spline_truc, aggc_taylor);
leg = legend('LSAP', 'Taylor', 'location', 'best');

figure('Name', 'Money Supply');
plot(spline_truc, m_lsap, spline_truc, m_taylor);
leg = legend('$m$ - LSAP', '$m$ - Taylor', 'location', 'best');

figure('Name', 'Inflation');
plot(spline_truc, pie_lsap, spline_truc, pie_taylor);
leg = legend('$\pi$ - LSAP', '$\pi$ - Taylor', 'location', 'best');

figure('Name', 'Output');
plot(spline_truc, y_lsap, spline_truc, y_taylor);
leg = legend('$y$ - LSAP','$y$ - Taylor', 'location', 'best');   